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(57) A method, apparatus, and article of manufacture for enhanced integration of a signal such as those 
generated by a chromatograph. In one embodiment an approximate baseline (1002) used for integration of a 
signal (101) is improved (1001) to reduce error. Alternate embodiments filter drift, noise, or both from a signal 
prior to integration or simulated distillation. The signal is processed according to a computer program stored 
in a memory. The methods include steps for subtracting (702) an approximate baseline (1002) from the signal 
(101), defining (703) a noise band (802-803) in the resulting difference (801), and forming a composite baseline 
(901) by substituting (700) signal data (101) for baseline data (1002) wherever the difference does not exceed 
the threshold (802-803), thus leaving the parts of the approximate baseline (1002) that correspond to signal 
peaks (102) unchanged (see eg. Fig 7). 
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SIGNAL ANALYSIS SYSTEM AND METHOD 

The present invention relates generally to signal analysis, especially two-dimensional 
signal analysis, and more particularly to a method, apparatus, and article of manufacture for 
the enhanced integration of time-varying signals which consist of one or more peaks rising 
above a background level containing noise and drift, such as those generated by 
chromatographic systems. 



A typical chromatographic analysis system includes an injection pon into which the 
sample is injected and mixed with an inert gas or liquid, a column through which the various 
dissolved components of the sample will travel at a rate related to the characteristics of the 
specific components, and a detector for measuring the concentration of each component as it 
exits the column and generating a time-varying signal corresponding to its concentration. The 
time between the injection of a sample and the detection of a specific component is called the 
retention time of that component. The time-varying signal is sent to a computing apparatus that 
integrates the signal by executing a computer program. The computing apparatus often provides 
a two-dimensional visual display of the signal, known in the art as a chromatogram. The x-axis 
of the chromatogram represents time, and the y-axis typically corresponds to the amplitude of the 
signal. 

A typical time-varying signal 101 generated by a gas chromatograph is shown in Figure 1. 
Peaks 102 in the signal correspond to the detection by the chromatograph of specific components 
of the chemical sample. The height of the peak above the background level 103 of the signal. 
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and the amount of area under the peak and above the background level, correlate to the amount of 
the respective component present in the sample. This information can be used for many 
applications, such as to check for the presence of, or to ensure the absence of, particular 
components in the sample, or to verify the proper concentration of the component in the sample. 

The computer program includes an integration method which detects peaks, determines 
the retention time 104 at which they occur, and determines their heights and areas. To determine 
peak height and peak area, the computer must first determine a baseline for each peak. The 
baseline is a reference level that corresponds to the background level of the signal when no peaks 
are present. The background level is typically a non-zero signal that may include noise (higher in 
frequency than the peaks) and drift (lower in frequency than the peaks; also known as wander). 

A chromatographer visually drawing the baseline under a particular peak would 
extrapolate it through the center of the noise on both sides of the peak with a slope corresponding 
to the drift. However, the presence of noise and drift in the background level make determining 
an accurate baseline difficult for a computer. As a result, prior art integration methods, in the 
presence of noise and drift, typically do not produce the same baseline as would a 
chromatographer looking at the signal visually. 

A first prior art integration method sets the baseline by looking at only a small region of 
the signal at a time. As shown in Figure 2, this method sets the baseline 201 too low where noise 
202 is present because the method finds the lowest points near both sides of the peak and draws 
the baseline 201 between these points. If, as shown in Figure 3. the detector produces a negative 
background signal disturbance 302 in the presence of two large peaks close together, this 
integration method will use the lowest point of the disturbance to set the baseline 301. Also, 
since the threshold level for signal detection has to be set conservatively in order to avoid 
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misidentifying background noise as peaks, some valid peaks may not be detected. This problem 
worsens as the signal-to-noise ratio decreases, as shown by the absence of any identified peaks 
for the signal 401 of Figure 4. 

A second prior art method extends the baseline back in time from the average background 
level that occurs towards the end of an analysis in order to compensate for the effects of venting 
the solvent peak. As shown in Figure 5. this method sets the baseline 501 too low in the region 
where background drift 504 occurs after the venting of the solvent peak 502. In addition, if the 
background continues to drift throughout the run. this method may misinterpret a hump in the 
drift as a false peak 505. 

A third prior art method establishes a baseline using curve fitting techniques. However, 
the complexity of curve fitting increases with the complexity of the baseline shape, which makes 
such a method difficult to implement and requires excessive computational power when the 
signal contains noise and drift. To reduce complexity, this method may alternatively divide the 
baseline into multiple segments to allow simpler, piecewise curve fitting. However, this 
modification introduces discontinuities between curve segments that affect the accuracy of 
integration results. 

A fourth prior art method filters noise by performing a fourier transform. However, it is 
difficult to filter out all the noise without distorting the peaks, since typically there is overlap in 
frequency content between the peaks and the noise. 

A chromatogram which displays an accurate baseline along with the signal is readily 
usable by an expert chromatographer. However, chromatograms are sometimes analyzed by non- 
chromatographers. for whom baseline drift and noise can cause confusion. A better display for 
non-experts would present a chromatogram from which the baseline noise and drift have been 



filtered, leaving only undistorted peaks and a flattened baseline. A filtered chromatogram could 
be accurately integrated using simple integration methods, and would provide more accurate 
results when analyzed using simulated distillation techniques. 

In the analysis of signals consisting of one or more peaks rising above a background level 
which contains noise, drift, or a combination of both, there has been a need for a simple method 
of determining an accurate baseline which eliminates noise and drift without distorting the peaks, 
in order to calculate accurate peak heights and areas. There has also been a corresponding need 
for filtering the signal to eliminate noise, drift, or both, prior to display, integration, or simulated 
distillation. 

The present invention may be implemented as a method of optimizing the baseline of a 
signal to account for noise and drift without introducing peak distortion. This method facilitates 
accurate measurements of peak height and peak area during integration, without undue 
complexity. The present invention may alternatively be implemented as a simple method of 
filtering the signal to attenuate noise, drift, or both. A chomatogram filtered according to this 
invention communicates information more clearly to non-chromatographers, can be accurately 
integrated by simple integration methods, and provides accurate results in simulated distillation 
analysis. The invention may be embodied in a program storage medium readable by a computer 
which contains program code that performs the baseline optimization and signal filtering 
methods. It may also be embodied in a computing apparatus having a computer program which 
performs the baseline optimization and signal filtering methods; a gas chromatographic apparatus 
may optionally be interfaced to the computing apparatus. 



A preferred first embodiment of the present invention is a method that produces an 
improved baseline from an approximate baseline by subtracting the approximate baseline from 
the signal, defming a noise band in the resulting difference, forming a composite baseline by 
substituting signal data for approximate baseline data at those segments of the baseline where the 
difference falls within the noise band (thus leaving the segments of the approximate baseline that 
correspond to the signal peaks unchanged), and smoothing the composite baseline to form the 
improved baseline. These steps may optionally be performed repetitively to further improve the 
baseline, using the improved baseline from the prior iteration as the approximate baseline for the 
subsequent iteration. A further option is to integrate the signal using the improved baseline to 
obtain the peak height and peak area. 

A preferred second embodiment of the present invention is a method that filters a signal 
to attenuate drift by smoothing the signal to form an approximate baseline, subtracting the 
approximate baseline from the signal, defining a noise band in the resulting difference, forming a 
composite baseline by substituting signal data for approximate baseline data at those segments of 
the baseline where the difference falls within the noise band (thus leaving the segments of the 
approximate baseline that correspond to the signal peaks unchanged), smoothing the composite 
baseline, and subtracting the smoothed composite baseline from the initial signal to form the 
filtered signal with drift attenuated. These steps may optionally be performed repetitively to 
further filter the drift, using the smoothed composite baseline from the prior iteration as the 
approximate baselme for the subsequent iteration. A further option is to integrate the filtered 
signal to obtain the peak height and peak area once the baseline has been improved using a 
simple integration algorithm. 
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A preferred third embodiment of the present invention is a method that filters a signal to 
, attenuate both noise and drift by smoothing the signal to form an approximate baseline, 
subtracting the approximate baseline from the signal, defining a noise band in the resulting 
difference, forming a composite baseline by substituting signal data for approximate baseline 
data at those segments of the baseline where the difference falls within the noise band (thus 
leaving the parts of the baseline that coaespond to the signal peaks unchanged), and subtracting 
the composite baseline from the initial signal to form the filtered signal with noise and drift 
attenuated. These steps may optionally be performed repetitively to further filter the noise and 
drift, by smoothing the composite baseline from the prior iteration and then using it as the 
approximate baseline for the subsequent iteration. A further option is to integrate the filtered 
signal to obtain the peak height and peak area once the baseline has been improved using a 
simple integration algorithm. Alternatively, a simulated distillation analysis may be performed 
on the filtered signal. 

A preferred fourth embodiment of the present invention is a method that filters a signal to 
attenuate noise by smoothing the signal to form an approximate baseline, subtracting the 
approximate baseline from the signal, defining a noise band in the resulting difference, forming 
an unsmoothed composite baseline by substituting signal data for approximate baseline data at 
those segments of the baseline where the difference falls within the noise band (thus leaving the 
parts of the baseline that correspond to the signal peaks unchanged), smoothing the composite 
baseline to form a smoothed composite baseline, and subu-acting the unsmoothed composite 
baseline from the initial signal while adding the smoothed composite baseline to the initial signal 
to form the filtered signal with noise attenuated. These steps may optionally be performed 
repetitively to further filter the noise, using the smoothed composite baseline from the prior 
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iteration as the approximate baseline for the subsequent iteration. 

Other aspects and advantages of the present invention will become apparent from the 
following detailed description, taken in conjunction with the accompanying drawings, illustrating 
by way of example the principles of the invention. The claims alone, not the preceding summary 
or the following detailed description, define the invention. 

Fig. 1 is an exemplary chromatogram showing a chromatographic signal including peaks 
and noise, and an approximate baseline. 

Fig. 2 is an exemplary chromatogram showing a chromatographic signal having 
background noise and an inaccurate baseline resulting from use of a prior art integration method. 

Fig. 3 is an exemplary chromatogram showing a chromatographic signal having a 
negative background disturbance and an inaccurate baseline resulting from use of a prior an 
integration method. 

Fig. 4 is an exemplary chromatogram showing a chromatographic signal with a low 
signal-io-noise ratio where no peaks are delected by a prior an integration method. 

Fig. 5 is an exemplary chromatogram showing a chromatographic signal including peaks, 
noise, and drift. The chromatogram shows an inaccurate baseline and a false peak resulting from 
use of a prior art integration method. 

Fig. 6 is a block diagram of a gas chromatographic instrumentation system suitable for 
the practice of the present invention. 

Fig. 7 is a simplified representation of an algorithm for performing baseline optimization 
in accordance with the present invention. 
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Fig. 8 is an exemplary chromatogram showing the baseline-adjusted signal generated by 
subtracting the approximate baseline from the signal of Fig. 1, and showing the threshold values 
defining the noise band for the signal. 

Fig. 9 is an exemplary chromatogram showing the composite baseline (before smoothing) 
generated by the first iteration of the baseline optimization algorithm operating on the signal of 
Fig. 1. 

Fig. 10 is an exemplary chromatogram showing, on an expanded scale relative to Fig. 1, 
the initial approximate baseline, the first pass improved baseline, and the second pass improved 
baseline generated by the baseline optimization algorithm operating on the signal of Fig. 1. 

Fig. 1 1 is an exemplary chromatogram showing the baseline-adjusted signal and the 
threshold values defining the noise band for the signal of Fig. 1, after the second iteration of the 
baseline optimization algorithm. 

Fig. 12 is a simplified representation of an algorithm for performing signal filtering to 
attenuate drift in accordance with the present invention. 

Fig. 13 is an exemplary chromatogram showing the filtered signal resulting from applying 
the signal filtering method of Fig. 12 to the signal of Fig. 5. 

Fig. 14 is a simplified representation of an algorithm for performing signal filtering to 
attenuate drift and noise in accordance with the present invention. 

Fig. 15 is an exemplary chromatogram showing the composite baseline (offset slightly 
downward for viewability) generated by applying the signal filtering method of Fig. 14 to the 
signal of Fig. 5. 

Fig. 16 is an exemplary chromatogram showing the filtered signal resulting from applying 
the signal filtering method of Fig. 14 to the signal of Fig. 5. 
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Fig. 17 is a simplified representation of an algorithm for performing signal filtering to 
attenuate noise in accordance with the present invention. 

Fig. 18 is an exemplary chromatogram showing the filtered signal resulting from applying 
the signal filtering method of Fig. 17 to the signal of Fig. 5. 

As shown in the drawings for purposes of illustration and as described in the following 
detailed description, the method, apparatus, and article of manufacture of the present invention 
improve the analysis of electronic signals consisting of one or more peaks greater in amplitude 
than a background level containing noise higher in frequency than the peaks, drift lower in 
frequency that the peaks, or both noise and drift. According to the present invention, integration 
and simulated distillation of these signals produce more accurate measurements, and the 
chromatograms of the analyzed signals can be understood more readily by non-experts. Prior 
methods of analyzing these signals have been unduly complex or have yielded less accurate 
measurements. 

Electronic signals usable with the present invention can be generated by a variety of 
instrumentation, and may represent a variety of physical, electrical, or chemical properties. 
Examples of such instrumentation include, but are not limited to, gas and liquid chromatographs, 
spectrum analyzers, oscilloscopes, and seismographs. Signals from some instruments may be 
generally bipolar, having sinusoidally-shaped peaks centered around the background level. The 
signals from chromatographs are generally unipolar, having Gaussian-shaped peaks rising above 
the background level (however, signals dropping below the background level, often called 
negative background disturbances, may occur as well). A gas chromatograph generates the 



signals according to the preferred practice of the invention, and therefore the detailed description 
of the invention will be directed to a gas chromatographic analytical system. However, it should 
be understood that the teachings herein are applicable to any electronic signals having the 
amplitude and frequency properties described above. 

Gas chromatographs usable in the practice of the present invention include the Hewlett- 
Packard 6890 series and 5890 series, while liquid chromatographs include the Hewlett-Packard 
1 100 series and 1050 series. 

System, Apparatus, and Program Storage Medium 
The preferred embodiment of an analytical instrument system 600 illustrated in Figure 6 
includes a gas chromatograph 610, and a computer system 620 having a signal input 621, 
processor 622. memory 623, interface 624, input device 627. and output device 628. In order to 
perform a chromatographic separation of a sample 601, a quantity of the sample 601 is injected 
into a fluid stream, preferably in the form of a pressurized carrier gas, via an injection port 61 1 
supplied with carrier gas stream 613. The injection port 61 1 provides a portion of the 
sample/carrier gas mixture 612 to a separation column 614. The column 614 is positioned within 
a temperature-controlled thermal chamber or oven 615. The carrier gas/sample combination 612 
passing through the column 614 is exposed to a temperature profile resulting from the operation 
of the heater 616 within the oven 614. 

As the carrier gas (containing the separated components of the sample) exits the column 
614. the presence of one or more of the constituent components of the sample is detected by a 
detector 617. The detector 617 can be any of the GC detectors known in the art, as long as it is 
capable of determining at least one physicochemical property of the carrier gas which exits the 
column 614. 
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Such GC detectors include the name ionization detector (FID), photoionization detector 
(PID). nitrogen phosphorous detector (NPD). flame photometric detector (FPD). thermal 
conductivity detector (TCD), atomic emission detector (AED), electrolytic conductivity detector 
(ELCD). and electron capture detector (ECD). Mass spectral detectors and infrared spectral 
detectors are also known. 

The detector output signal 618 may be either in analog or digital electronic form. The 
detector output signal is then received by the signal input 621 for communication to the processor 
622. If the signal is in analog form, the signal input may include an analog-to-digital converter 
(not shown) to convert the signal to a computer-readable form. The converted detector output 
signal is typically stored in the memory 623. Preferably, the detector output signal is provided in 
the form of data representative of one or more peaks each corresponding to detection of a 
constituent component of the sample, the peak(s) rising above a background signal level 
containing noise, drift, or both noise and drift. The operations performed by the computer system 
620 on the signal data in accordance with the present invention will be described subsequently. 

The processor 622 may include computing devices amenable to the practice of this 
invention, e.g., one or more computing devices such as computers, microprocessors, 
microcontrollers, switches, logic gates, or any equivalent logic device capable of performing the 
computations described herein below. The processor 622 preferably is coupled with the interface 
624. input device 627, and output device 628. Input devices preferably include one or more of a 
keyboard, keypad, pointing device, or remote computer (not shown i for inputting operating 
condition parameters, system calibration data, and the like. Output devices preferably include 
alphanumeric or video displays, printers, or a remote computer (not shown). 
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The preferred compuier system 620 further includes memory 623 in the form of volatile 
and non-volatile devices in which computer programs 626, input and output information, 
operating condition parameters, and system information can be stored and retrieved. 
Contemplated computer programs 626 include method embodiments of the present invention. 
Contemplated memory devices include, but are not limited to, ROM, RAM, floppy disk, hard 
disk, and CD-ROM. Operating commands, detector response attributes, column temperature 
programs, and other information necessary to perform a chromatographic analysis may be entered 
into the processor 622 by way of an input device 627 or retrieved from memory 623. 

The chromatogram, peak heights, and peak areas calculated by the processor 622 for the 
signal may be transmitted to an output device 628 and displayed or printed. The interface 624 
may further comprise network and bus system (input/output or I/O) controllers, isolation devices, 
clocks, and other related electronic components for performing control, processing, and 
communication tasks other than those described herein. A remote computing apparatus (not 
shown) may connect to the interface 624 and cooperate with the computer system 620 to 
implement the present invention as a client/server application. 

The method embodiments of the present invention are operable in the system 600, and 
may be contained in the memory 623. 

Signal Data Representation and Storage 
The signal data, herein referred to mathematically as s(x), is stored in the memory 623 in 
a format usable for two-dimensional data. The dependent variable s(.x) represents the signal in a 
chromatogram. and the independent variable x represents time. In a first embodiment where the 
values of x represent a fixed interval of time, s(x) may be stored in a data stnicture comprising a 
one-dimensional array of signal data points. In a second embodiment where the values of x 
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represent a variable interval of time, a data structure comprising an array of ordered pairs (x,s(x)) 
may be used. Data structures of the corresponding type may be used to contain the other 
dependent variable data disclosed in the remainder of the detailed description, such as baselines 
b(x) and baseline-adjusted signals a(x). 

Smoothing Algorithm 

Several embodiments of the present invention include a method step for smoothing a 
signal or a baseline. The purpose of this step is to attenuate high frequencies (noise and peaks) in 
the signal or baseline while having little or no effect on low frequencies (drift). Smoothing 
conserves area; when a peak gets smoothed, its area gets redistributed over a wider range. 
Visually, the peak becomes lower in amplitude but extends over a wider period of time. 

A variety of high-frequency filtering algorithms may be utilized to perform the smoothing 
step. The general approach to smoothing is to replace each of the original data points with the 
weighted average of it and its nearest neighbors. The number of points averaged, and the 
weighting factors applied to each point, determine the frequency-response curve. There is one 
complication to any smoothing algorithm: data points near the start and end of the data do not 
have the required number of nearest neighbors on one side, so their smoothed values cannot be 
calculated in the usual way. Often the solution to this problem is to have fewer points in the 
smoothed data than were in the original data. The first and last (n-l)/2 points (where n is the 
number of points averaged) are simply missing. Where multi-pass smoothing is applied, more 
data will being lost at each subsequent pass. 

The preferred smoothing algorithm uses a simple "boxcar" filter in which all of the 
weighting factors are the same and the only variable is the number of points averaged. To 
mitigate the data distortion caused by boxcar filtering, a second pass is performed with fewer 
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points averaged. This achieves high-frequency rejection with less distortion, and is functionally 
equivalent to smoothing with unequal weighting factors. To avoid data loss at the start and end 
of the signal or baseline, the algorithm extrapolates additional data points prior to smoothing 
which are included in the filtering. 

An alternate smoothing embodiment having unequal weighting factors uses aSavitzky- 
Golay (least-squares) filter. The variables in least-squares smoothing are the number of points 
averaged and the order of equation to which the data is to be fit. Although least-squares 
smoothing is more complex than boxcar smoothing, it can achieve a much better tradeoff 
between high-frequency attenuation and data distortion in a single pass. Data loss at the start and 
end of the signal or baseline can be avoided by using weighting factors that calculate the least- 
squares value at a point other than the center point for each data point near the ends of the data. 

Where an embodiment of the present invention calls for smoothing a signal that contains 
large, broad peaks to form an initial or approximate baseline, a more accurate initial baseline may 
be constructed after smoothing by forming a composite consisting of smoothed signal at those 
values of x where the signal approximately represents the background, and straight line segments 
between the endpoints of the smoothed signal to bridge the gaps where the signal approximately 
represents peaks. In order to determine where the signal represents peaks, the steps of forming a 
baseline-adjusted signal 702 and measuring the noise to set threshold values 703 may be 
performed, as shown in Figure 7 and described subsequently. 

Noise Measuring Algorithm 

Several embodiments of the present invention include a method step for measuring the 
noise in a baseline-adjusted signal to define positive and negative threshold values, tp and In. 
The general approach to noise measurement involves dividing the signal up into smaller 
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contiguous groups of data points and measuring the noise in each segment. In each segment, 
noise can be measured as either peak-to-peak or RMS noise. Peak-to-peak noise is the difference 
between the largest and the smallest values within the segment. RMS is the standard deviation of 
all of the data values within the group. Peak-to-peak is faster and easier to measure, but RMS is 
more accurate. Peak-to-peak noise is roughly six times larger than RMS noise for the same data. 
The preferred embodiment measures RMS noise. The thresholds ip and ts are set to plus and 
minus a factor times the RMS noise respectively. The preferred embodiment uses a factor of 2.0. 

In measuring noise, the algorithm must ensure that it does not confuse peaks or drift with 
noise. In the preferred embodiment of noise measurement, any drift in the input signal has 
previously been subtracted from the baseline-adjusted signal in which the noise is measured. In 
order to exclude peaks, the preferred embodiment finds those segments with the smallest non- 
zero noise (background noise is assumed to be much lower than peaks). Once the noise 
measurements have been made for all segments, the smallest value can be chosen as the baseline 
noise. Since a particular detector may saturate above a certain signal level, causing it to have 
erroneously low noise, the preferred embodiment rejects any measured noise value of zero and 
selects the smallest non-zero value as the noise. 

A funher aspect of the noise measurement algorithm involves which groups of points are 
used to perform the measurement. The grouping can be as simple as dividing up the data points 
into some number of approximately equal groups (for example, one hundred groups of n/ 100 
each, where n is the number of points in the signal) and calculating the noise of each group. A 
more sophisticated approach measures noise on more than 100 groups but still includes 1/100 of 
the data points in each group, for example the first group might be points 1 to 100, the second 
group points 51 to 150, etc. The preferred embodiment uses groups of data at one -data-point 
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increments of the starting data point (for example, points I to 100. 2 to 101, 3 to 102, etc.). 
Baseline Optimization Method 
A preferred first embodiment of the present invention is a novel method of optimizing a 
baseline used for integration of a signal. The general approach is to begin with an initiaJ 
approximate baseline that approximately tracks the background drift of the signal but has an error 
component, and improve the accuracy of this baseline by forming a composite baseline which 
tracks the signal noise and then smoothing the composite baseline to center the baseline within 
the noise. The error component of this improved baseline is smaller than that of the initial 
baseline. 

As illustrated by applying the method of Figure 7 to the signal data of Figure 1 by way of 
example, the first embodiment begins with a initial step 701 in which an approximate baseline 
b(x) 105 is generated. The approximate baseline may be generated by smoothing the signal s(x), 
or by other methods which produce a baseline that generally tracks the background drift of the 
signal. The approximate baseline contains inaccuracies, particularly in the region of signal peaks 
and negative background disturbances, due to the area conservation of the smoothing algorithm. 

Once the approximate baseline has been obtained, the following optimization steps can be 
performed one or more times to optimize the baseline. The first step 702 in optimization is to 
subtract the baseline from the signal, the difference forming a baseline-adjusted signal a(x) 801 
centered around zero counts, as shown in Figure 8. The baseline-adjusted signal shows the 
deviation of the signal from the baseline; for example, at a value of x where the basehne b(x) 
equals the signal s(x), the baseline-adjusted signal a(x) equals zero counts. Since the smoothed 
baseline from the prior step tracked the drift, drift will be attenuated in the baseline-adjusted 
signal. 
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The second optimizing step 703 is to calculate the noise in the baseline-adjusted signal 
a(x) 801 and set positive (tp) 802 and negative (In) 803 threshold values outside the range of, and 
surrounding, the noise. Signal data s(x) at those values of x where the baseline-adjusted signal 
falls within the noise band defined by the thresholds approximately represent signal background. 
Signal data at those values of x where the baseline-adjusted signal falls outside this band 
approximately represent signal peaks or negative background disturbances. 

The third step 704 is to substitute data from the signal s(x) into the baseline b(x) at those 
values of x where the baseline-adjusted signal a(x) falls within the noise band defined by the 
threshold values. After this step has been performed, the baseline 901 in Figure 9 becomes a 
composite of signal data 902 at those values of x where the signal approximately consists of 
background, and previous baseline 903 at those values of x where the signal approximately 
consists of peaks or negative background disturbances. As will be appreciated by comparing 
Figure 9 with Figure 1, the composite baseline 901 approximates signal 101 data, with peaks 102 
and disturbances 106 outside the noise band clipped. Replacing baseline data b(x) with signal 
data s(x) has the effect of removing from the composite baseline 901 much of the area of the 
peaks which was included in the baseline 105 following smoothing 701. 

The fourth step 705 is to smooth the composite baseline 901 to form the improved 
baseline 1001 which has a reduced error component. Since signal data representing peaks and 
disturbances outside the threshold values have been removed from the composite baseline 901, 
they do not distort the first pass improved baseline 1001 in Figure 10 resulting from the 
smoothing step 705. It will be appreciated that the first pass improved baseline 1001 is closer to 
a baseline centered through the noise and drift than the initial baseline 1002, and is closer to a 
baseline that a chromatographer looking at the chromatogram would draw. 
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The preceding optimizing steps 702. 703, 705. and 705 may be performed iteratively 706 
to further improve the baseline, substituting the improved baseline generated by the previous 
iteration as the approximate baseline for the next iteration. Because the first pass improved 
baseline 1001 is more accurate than the initial baseline 1002. the deviations from zero counts in 
the second pass baseline-adjusted signal 1 101 of Figure 1 1 are further reduced, allowing second- 
pass thresholds tp 1 102 and tN II 03 to be set tighter relative to first-pass thresholds tp 802 and In 
803, which in turn results in a more accurate second-pass composite baseline. The second 
iteration yields the second pass improved baseline 1003 of Figure 10. For a fairly simple signal 
with small peaks and disturbances, good results are possible after only one or two iterations, but 
to handle the more general case the preferred embodiment iterates eight times. Additional 
iterations of the optimization steps yield further improvement, at the cost of increased processing 



time. 



After the baseline has been optimized, the signal may be integrated via a prior art 
integration algorithm that uses this optimized baseline to define the time and signal level to use 
for the starting and ending baseline points for each peak so as to more accurately calculate peak 
heights and areas. 

Signal Filtering Method to Attenuate Drift 

As illustrated by applying the method of Figure 12 to the signal data of Figure 5 by way 
of example, a preferred second embodiment of the present invention is a novel method of 
filtering a signal 503 to attenuate background drift 504. The general approach used in the second 
embodiment is to generate a smoothed composite baseline which tracks the signal drift, and then 
subtract the baseline from the signal to attenuate the drift. The similarities between this method 
and the baseline optimization method will be appreciated by those skilled in the art. 
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This embodiment begins with a step 1201 in which an initial baseline bCx) is generated by 
smoothing the signal s(x) 503. The smoothing operation generates a baseline that approximately 
tracks low-frequency phenomena such as background drift but not higher-frequency phenomena 
such as background noise and peaks. 

The second step 1202 is to subtract the baseline b(x) from the signal s(x) 503, the 
difference forming a baseline-adjusted signal a(x) centered around zero counts. Since the 
smoothed baseline b(x) tracks drift in the signal but not peaks or noise, it will be appreciated that 
drift in the baseline-adjusted signal a(x) is attenuated while peaks and noise remain. 

The third step 1203 is to calculate the noise in the baseline-adjusted signal a(x) and set 
positive (tp) and negative (ts) threshold values outside the range of. and surrounding, the noise. 
Signal data at those values of x where the baseline-adjusted signal falls within the noise band 
defined by the thresholds approximately represent background. Signal data at those values of x 
where the baseline-adjusted signal falls outside this band approximately represent peaks. 

The fourth step 1204 is to substitute data from the signal s(x) into the baseline b(x) at 
values of x where the baseline-adjusted signal a(x) falls within the noise band defined by the 
threshold values. Following the performance of this step, the baseline becomes a composite of 
signal data at those values of x where the signal approximately consists of background, and 
previous baseline at those values of x where the signal approximately consists of peaks. The 
composite baseline approximates signal data with any peaks outside the noise band clipped. It 
will be appreciated that substituting signal data into the baseline makes the composite baseline 
approximately track the drift and noise in the signal. 



20 

The fifth step 1205 is to smooth the composite baseline b(x). The smoothing operation 
results in a baseline b(x) that still approximately tracks background drift, but the effects of noise 
get smoothed out and attenuated. 

The second, third, founh. and fifth steps 1202, 1203, 1204, and 1205 above may be 
performed iteratively 1207 to improve how well the baseline b(x) tracks the background drift, by 
using the smoothed composite baseline resulting from the previous iteration as the approximate 
baseline for the next iteration. For a fairly simple signal with small peaks and disturbances, good 
results are possible after only one or two iterations, but to handle the more general case the 
preferred embodiment iterates eight times. Additional iterations further attenuate drift, at the 
cost of increased processing time. 

After the final iteration has been performed, the last step 1206 is to subtract the baseline 
b(x) from the signal s(x) to form the filtered signal 1301 shown in Figure 13. It will be 
appreciated that this results in a signal 1301, centered around zero counts, in which background 
drift has been attenuated but background noise has not. The heights of the peaks 1302 have been 
adjusted to compensate for drift, as can be seen by comparison to the peaks 506. Peak area is 
conserved. 

The filtered signal may further be displayed as a chromatogram, or it may be integrated 
via a prior art integration algorithm to calculate peak heights and areas. Since drift has been 

filtered from the signal, integration is simplified in that a signal level of zero counts may be used 
for the starting and ending baseline points for all peaks. In addition, false peaks due to drift are 
elinninated. as can be seen by the absense in Figure 13 of the false peak 505. 
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Signal Filtering Method to Attenuate Noise and Drift 
As illustrated by applying the method of Figure 14 to the signal data of Figure 5 by way 
of example, a preferred third embodiment of the present invention is a novel method of filtering a 
signal 503 to attenuate both noise 507 and drift 504. The general approach is to generate a 
smoothed baseline which tracks the signal drift, then create a composite baseline which tracks the 
signal noise, and finally subtract the composite baseline from the signal to attenuate the drift and 
noise. The similarities between this method and the signal filtering method to attenuate drift 
described previously will be appreciated by those skilled in the an. 

This embodiment begins with a step 1401 in which an initial approximate baseline b(x) is 
generated by smoothing the signal s(x) 503. The smoothing operation generates a baseline that 
tracks low-frequency phenomena such as drift but not higher-frequency phenomena such as noise 
and peaks. 

The second step 1402 is to subtract the baseline b(x) from the signal s(x) 503, the 
difference forming a baseline-adjusted signal a(x) centered around zero counts. Since the 
smoothed baseline b(x) tracks drift in the signal but not peaks or noise, it will be appreciated that 
drift in the baseline-adjusted signal a(x) is attenuated while peaks and noise remain. 

The third step 1403 is to calculate the noise in the baseline-adjusted signal a(x) and set 
positive (tp) and negative (In) threshold values outside the range of. and surrounding, the noise. 
Signal data at those values of x where the baseline-adjusted signal falls within the noise band 
defined by the thresholds approximately represent background. Signal data at those values of x 
where the baseline-adjusted signal falls outside this band approximately represent peaks. 

The fourth step 1404 is to substitute data from the signal s{x) into the baseline b(x) at 
values of X where the baseline-adjusted signal a{x) falls within the noise band defined by the 
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threshold values. Following the performance of this step, the baseline becomes a composite of 
signal data at those values of x where the signal approximately consists of background, and 
previous baseline at those values of x where the signal approximately consists of peaks. The 
composite baseline approximates signal data with any peaks outside the noise band clipped. It 
will be appreciated that substituting signal data into the baseline makes the composite baseline 
approximately track the drift and noise in the signal. 

The second, third, and fourth steps above may be performed iteratively 1406 to further 
improve how well the composite baseline tracks the background drift and noise. If the method is 
performed iteratively. the composite baseline b(x) generated by the fourth step 1404 of the prior 
iteration will be smoothed 1405 and then will replace the approximate baseline used for 
performing the subtraction step 1402 of the next iteration. For a fairly simple signal with small 
peaks and disturbances, good results are possible after only one or two iterations, but to handle 
the more general case the preferred embodiment iterates eight times. Additional iterations 
further attenuate drift and noise, at the cost of increased processing time. Figure 15 shows the 
composite baseline 1501 produced by the fmal iteration, overlaid with the signal 1502 (for 
purposes of clarity, the baseline is illustrated offset slightly downward). It will be appreciated 
that the baseline closely tracks the signal everywhere but where peaks are present. 

Once iteration has been completed and the final composite baseline b(x) has been 
obtained, the composite baseline is subtracted 1407 from the signal s(x). As shown in Figure 16, 
the background drift and noise are attenuated; in fact, after sufficient iterations the resulting 
signal s(x) 1601 will be zero counts everywhere except where peaks 1602 exist. The height of 
the peaks 1602 above the baseline is unchanged, as can be seen by comparison to the peaks 506, 
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and peak area is conserved. In addition, false peaks due to drift are eliminated, as can be seen by 
the absence in Figure 16 of the false peak 505. 

The drift-adjusted signal may further be displayed as a simple chromaiogram, or it may be 
integrated via a prior art integration algorithm to calculate peak heights and areas. Since drift 
and noise have been filtered from the signal, integration is simplified. A signal level of zero 
counts is used for the starting and ending baseline points of all peaks, and a simple peak detector 
will suffice since any positive deviation from zero counts indicates a peak. In addition, false 
peaks due to drift are eliminated, as can be seen by the absence in Figure 16 of the false peak 
505. 

Another application of this filtering method is in simulated distillation, which does not 
involve integration. In simulated distillation the signal is divided up into small segments called 
area slices. Each area slice corresponds to a certain boiling range of a hydrocarbon mixture. 
Which peaks make up the area within a slice is not important, since all peaks in the slice have 
approximately the same boiling point; rather, it is the total area under the peaks of each slice that 
are measured. Significant errors result from simulated distillation when noise and drift are 
significant relative to the signal. However, by using this filtering method to remove drift and 
noise from the signal prior to performing simulated distillation, the area slices will only contain 
the area of the peaks. This technique is only applicable where there are discrete peaks, as in 
lower-boiling mixtures. Higher boiling ranges tend to produce "humps" of highly merged peaks 
and the filtering method will treat the humps as drift and subtract out much of their area. 

Signal Filtering to Attenuate Noise 

As illustrated by applying the method of Figure 17 to the signal data of Figure 5 by way 
of example, a preferred fourth embodiment of the present invention is a novel method of filtering 



a signaJ 503 to attenuate noise 507. The general approach is to form a composite baseline which 
tracks noise and drift, then smooth the composite baseline to form an intermediate baseline 
which tracks only drift, subtract the composite baseline from the signal to attenuate both noise 
and drift, and finally add the intermediate baseline back in to reconstitute the drift. Those skilled 
in the art will appreciate that this method combines elements of the two signal filtering methods 
described above. 

This embodiment begins with a step 1701 in which an initial baseline b(x) is generated by 
smoothing the signal s(x) 503. The smoothing operation generates a baseline that tracks low- 
frequency phenomena such as drift but not higher-frequency phenomena such as noise and peaks. 

The second step 1702 is to subtract the baseline b(x) from the signal s(x), the difference 
forming a baseline-adjusted signal a(x) centered around zero counts. Since the smoothed 
baseline b(x) tracks drift in the signal but not peaks or noise, it will be appreciated that the drift 
in the baseline-adjusted signal a(x) is attenuated while peaks and noise remain. 

The third step 1703 is to calculate the noise in the baseline-adjusted signal a(x) and set 
positive (tp) and negative (ts) threshold values outside the range of, and surrounding, the noise. 
Signal data at those values of x where the baseline-adjusted signal falls within the noise band 
defined by the thresholds approximately represent background. Signal data at those values of x 
where the baseline-adjusted signal falls outside this band approximately represent peaks. 

The founh step 1704 is to substitute data from the signal s(x) into the baseline b(x) at 
values of X where the baseline-adjusted signal a(x) falls within the noise band defined by the 
threshold values. Following the performance of this step, the baseline becomes a composite of 
signal data at those values of x where the signal approximately consists of background, and 
previous baseline at those values of x where the signal approximately consists of peaks. The 
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composite baseline approximates signal data with any peaks outside the noise band clipped. It 
will be appreciated that substituting signal data into the baseline makes the composite baseline 
approximately track the background drift and noise in the signal. 

The fifth step 1705 is to smooth the composite baseline b(x) to form an intermediate 
baseline b2(x) that still approximately tracks background drift, but no longer tracks background 
noise due to effects of the smoothing operation. 

The second, third, fourth, and fifth steps above may be performed iteratively 1706 to 
further improve how well the composite baseline b(x) and the intermediate baseline b2(x) track 
the background. If the method is performed iteratively. the intermediate baseline b2(x) generated 
by the fifth step 1705 of the prior iteration will replace the baseline b(x) 1707 for the subtraction 
step 1702 of the next iteration. For a fairly simple signal with small peaks and disturbances, 
good results are possible after only one or two iterations, but to handle the more general case the 
preferred embodiment iterates eight limes. Additional iterations further attenuate noise, at the 
cost of increased processing time. 

Once iteration has been completed and the final composite baseline b(x) and intermediate 
baseline b2(x) have been obtained, the composite baseline is subtracted 1708 from the signal s(x) 
to attenuate the background drift and noise, and the intermediate baseline b2(x) is added 1708 to 
the signal s(x) to add the background drift back in. As shown in Figure 18. the background noise 
is attenuated while the drift 1803 is generally unaffected. As a result, the amplitude of peaks 
1802 is substantially the same as the peaks 506 in the original signal 503. It will be appreciated 
that because of the addition of the background drift, the resulting signal s(x) 1801 will not be 
centered around zero counts, but around the level of the original background drift 1803. 
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The drift-adjusted signal may further be displayed as a chromatogram which shows the 
drift in the signal without the masking effects of noise. The signal may also be integrated via a 
prior art integration algorithm to calculate peak heights and areas. 

Conclusion 

From the foregoing it will be apparent that the invention provides novel and 
advantageous methods of optimizing the baseline of a signal to properly account for the 
effects of noise and drift, and of filtering drift, noise, or both from a signal prior to 
performing subsequent processing operations. An accurate baseline ensures proper peak 
recognition and produces accurate integration results for peak height and peak area. Pre- 
filtering the signal to remove drift and noise provides a visually cleaner display, allows 
integration to be performed by simpler and faster algorithms, and allows accurate 
simulated distillation to be performed. 

Although several specific embodiments of the invention have been described and 
illustrated in the area of gas chromatography, the invention is not to be limited to the 
specific forms, methods, or arrangements of parts so described and illustrated. For 
example, the embodiments of the invention described herein may operate upon signals 
which are generated by other electronic instruments such as spectrum analyzers and 
seismographs. The invention is defined by the claims. 
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aAIMS 

1 . A method for analyzing an input signal having a background level, comprising the method 
steps of: 

acquiring the input signal (101); 

subtracting (702) an approximate baseline (1002) having a first error component from the 
input signal (101) to forai a baseline-adjusted signal (801) having baseline-adjusted noise; 

calculating (703) first (802) and second (803) threshold values surrounding the range of 
the baseline-adjusted noise and definite a noise band between the first and second threshold 
values; and 

substituting (704) the input signal (101) for the approximate baseline (1002) wherever the 
baseline-adjusted signal (801) falls within the noise band to form a composite baseline (901). 

2. The method for analyzing an input signal as claimed in claim 1, fiirther comprising the method 
step of: 

smoothing (701) the input signal (101) to form the approximate baseline (1002). 

3. The method for analyzing an input signal (503) as claimed in claim 1, wherein the background 
level includes a first drift component, fiirther comprising the method steps of: 

smoothing (1205) the composite baseline to form an improved baseline; and 
subtracting (1206) the improved baseline from the input signal to form a filtered signal 
(1301) having a second drift component attenuated relative to the first drift component. 

4. The method for analyzing an input signal as claimed in claim 3, further comprising the 
method step of replacing the approximate baseline with the improved baseline, and wherein the 
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method steps of subtracting the approximate baseline, calculating first and second threshold 
values, substituting the input signal into the approximate baseline, smoothing the composite 
baseline, and replacing the approximate baseline are perfonned iteratively (1207), whereby the 
second drift component is further attenuated with each iteration. 

5. The method for analyzing an input signal (503) as claimed in claim 1. wherein the background 
level includes a first drift component and a first noise component, further comprising the method 
step of subtracting (1407) the composite baseline fi-om the input signal to form a filtered signal 
(1601) having a second drift component attenuated relative to the fiist drift component, and a 
second noise component attenuated relative to the first noise component 

6. The method for analyzing an input signal as claimed in claim 5, fiirther comprising the 
method steps of: 

smoothing (1405) the composite baseline to form an improved baseline, and 

replacing the approximate baseline with the improved baseline; 

wherein the method steps of subtracting the approximate baseline, calculating first and 
second threshold values, substituting the input signal into the a^Hoximate baseline, smoothing 
the composite baseline, and replacing the approximate baseline are perfonned iteratively (1406), 
whereby the second drift component and second noise component are fiirther attenuated with 
each iteration. 
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7. The method for analyzing an input signal as claimed in claim 5, wherein the input signal is 
chromatographic data, further comprising the method step of performing simulated distillation 

calculations on the filtered signal. 

8. The method for analyzing an input signal (503) as claimed in claim 1, wherein the background 
level includes a first noise component, further comprising the method steps of: 

smoothing (1705) the composite baseline to fonn an intermediate baseline; and 
forming (1 708) a filtered signal (1 801) having a second noise component attenuated 

relative to the first noise component by subtracting the composite basdine from, and adding the 

intermediate baseline to, the input signal. 

9. The method for analyzing an input signal as claimed in claim 8, findier comprising die 
method step of replacing (1707) the approximate baseline with the intennediate baseline, and 
wherein the method st^ of subtracting the approximate baseline, calculating first and second 
threshold values, substituting the input signal into the approximate baseline, smoothing the 
composite baseline, and replacing the approximate baseline are performed iteratively (1706), 
whereby the second noise component is further attenuated with each iteration. 

1 0. The method for analyzing an input signal as claimed in claim 1 , further comprising the 

method step of: 

smoothing (705) the composite baseline (901) to form an improved baseline (1001) 
having a second error component smaller than the first error component 



1 1 . The method for analyzing an input signal as claimed in claim 10. further comprising the 
method step of replacing the approximate baseline (1002) with the improved baseline (IQOl), 
wherein the method steps of subtracting (702) the approximate baseline (1002), calculating (702) 
first (802) and second (803) threshold values, substituting (704) the input signal (101) into the 
approximate baseline (1002), smoothing (705) the composite baseline (901). and replacing the 
approximate baseline (1002) are perfonned iteratively (706), and when:by the second error 
component is further attenuated with each iteratioa 

12. The method for analyzing an input signal as claimed in claim 1. ^in the input signal is 
chromatographic data. 

13. The method for analyzing an input signal as claimed in claim 10, further comprising the 
method step of integrating (707) the input signal (101) using the improved baseline (1001). 

14. The method for analyzing an input signal as claimed in claim 10. wherein the approximate 
baseline (1002) further comprises a data structure b(x), the input signal (101) further comprises a 
data structure s(x), the baseline-adjusted signal (801) further comprises a data structure a(x). and 
wherein the method step of substituting (704) the input signal (101) for the approximate baseline 
( 1 002) further comprises substituting (704) the value of s(x) into b(x) at aU values of x where 
a(x) falls within the noise band. 

1 5. A program storage medium (623) readable by a computer (622), tangibly embodying a 
program (626) of instructions executable by the computer to perform method steps for improving 
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the accuracy of an approximate baseline (1002) for an input signal (101), the approximate 
baseline (1002) having a first error component, and the input signal having a background level, 
comprising the method steps of: 

subtracting (702) the approximate baseline (1002) from the input signal (101) to form a 
baseline-adjusted signal (801) having baseline-adjusted noise; 

calculating (703) first (802) and second (803) threshold values surrounding the range of 
the baseline-adjusted noise and defining a noise band between the first and second threshold 
values; 

substituting (704) the input signal (101) for the ^proximate baseline (1002) wherever the 
baseline-adjusted signal (801) falls within (902) the noise band to form a composite baseline 
(901); and 

smoothing (705) the composite baseline (901) to form an improved baseline (1001) 
having a second error component smaller than the first error component. 

16. A program storage medium (623) readable by a computer (622), tangibly embodying a 
program (626) of instructions executable by the computer to perform method steps for filtering 
an input signal (503) having a background level, comprising the method steps of: 

smoothing (1201) the input signal (503) to form an approximate baseline; 

subtracting (1202) the approximate baseline from the input signal (503) to form a 
baseline-adjusted signal having baseline-adjusted noise; 

calculating (1203) first and second threshold values surrounding the range of the 
baseline-adjusted noise and defining a noise band between the fiirst and second threshold values; 
and 



subsututing (1204) the input signal (503) for the approximate baseline wherever the 
baseline-adjusted signal falls within the noise band to form a composite baseline. 

17. The program storage medium (623) readable by a computer (622) as claimed in claim 16, 
wherein the background level includes a first drift component, fimher comprising the method 
steps of: 

smoothing (1205) the composite baseline to form an improved baseline; and 
subtracting (1206) the improved baseline bom the input signal (503) to form a filtered 
signal (1301) having a second drift component attenuated relative to the first drift component. 

18. Tlje program storage medium (623) readable by a computer (622) as claimed in claim 16, 
wherein the background level includes a first drift component and a first noise component, 
fiirther comprising the method step of subtracting (1407) the composite baseline from the input 
signal (503) to fcmi a filtered signal (1601) having a second drift component attenuated relative 
to the first drift component, and a second noise component attenuated relative to the first noise 
component. 

19. The program storage medium (623) readable by a computer (622) as claimed in claim 16. 
wherein the background level includes a first noise component, fiirther comprising the method 
steps of; 

smoothing (1705) the composite baseline to form an intermediate baseline; and 
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fonning (1708) a filtered signal (1801) having a second noise component attenuated 
relative to the first noise component by subtracting the composite baseline fi-om, and adding the 
intcnnediate baseline to, the input signal (503). 

20. A data analysis system, comprising: 

an input (621) for receiving an input signal (618) having at least one peak and fiirther 
having a background level including background noise; and 

a processor (622) in communication with the input (621) wiiich 

subtracts (702) an approximate baseline bom the input signal to form a baseline- 
adjusted signal having baseline-adjusted noise, 

calculates (703) first and second direshold values sxirrounding the range of the 
baseline-adjusted noise thus defining a noise band between the first and second threshold values, 
and 

substitutes (704) the input signal for the approximate baseline wherever the 
baseline-adjusted signal falls within the noise band to form a composite baseline. 

21 . The data analysis system of claim 20, fiorther comprising a chromatograph (6 1 0) connected 
to the input (621), the chromatograph (610) adapted for receiving a chemical sample (601) 
having one or more components, the chromatograph (610) fiirther having a column (614) for 
separating the chemical sample (601) into the components and a detector (617) responsive to the 
components for generating the input signal (618). 
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22. The data analysis system of claim 20, ftmher comprising an output device (628), wherein ti 

processor (622) further 

smoothes (705) the composite baseline (901) to form an improved baseline (lOOI), 
integrates (707) the input signal (618) using the improved baseline (1001), and 
transmits the integrated signal to the output device (628). 

23. A method of signal analysis substantially as herein described with reference to Figs 
6 to 18 of the accompanying drawings. 

24. A program storage medium substantially as herein described with reference to Figs 
6 to 18 of the accompanying drawings. 



25. A data analysis system substantially as herein described with reference to Figs 6 to 
18 of the accompanying drawings. 
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